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ABSTRACT 



Context. A massive objects such as a globular cluster passing through the disk of a galaxy can trigger star formation. 
Aims. We test the hypothesis that the most massive globular cluster in the Galaxy, u Centauri, which crossed the disk 
approximately 24 ± 2Myr ago, may have triggered the formation of the open clusters Stephenson 2 and BDSB 122. 
Methods. The orbits of tj Centauri, Stephenson 2 and BDSB 122 are computed for the three-component model of 
Johnston, Hernquist & Bolte, which considers the disk, spheroidal and halo gravitational potentials. 
Results. With the re-constructed orbit of a; Centauri, we show that the latest impact site is consistent, within important 
uncertainties, with the birth-site of the young massive open clusters BDSB 122 and Stephenson 2. Within uncertainties, 
this scenario is consistent with the time-scale of their backwards motion in the disk, shock wave propagation and delay 
for star formation. 

Conclusions. Together with open cluster formation associated to density waves in spiral arms, the present results are 
consistent with the idea that massive globular clusters as additional progenitors of open clusters, the massive ones in 
particular. 
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1. Introduction 

Disk-stability criteria and impact assumptions suggest that 
the passage of a globular cluster (GC) can trigger a bubble 
or wave of sel f-propagating star formation within the d isk of 
the Galaxy (IWalhn. Higdon fc Stavelev-SmithI IT996I ). The 
initial mechanical perturbation produces a local enhance- 
ment of the interstellar medium (ISM) density, from which 
localised star formation may occur. Subsequently, clus- 
tered star formation may happen along the border of a 
radially expanding density wave or io nisation front (e.g. 
Soria ct al. 200 5, - hereafte r SCP05; lElme green fc Ladal 
19771: iWhitworth et al.lll994[ ). The expanding bubble is ca- 
pable of compressing the neutral ISM above the stability 
criterion against gravitational collapse. Alternative star- 
formation triggering mechani sms are the infall of a high- 
velocity HI clou d (Elmcgrc en. Efremov fc Lar sen 200^ 
Larsen et al.l l2002f ). or hypernova explosions ijPaczvhskil 
19981) . 

Prominent, isolated star- forming bubbles have been ob- 
served in external galaxies. A bubbl e with » 600 pc i n 
diameter was detected in NGC 6946 (jLarsen et al.l l2002| ) , 
containing a young super star cluster and at least 12 
surrounding young clusters, the latter being compara- 
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ble to the most luminous Galactic OCs. The trigger- 
ing mechanism in NGC 6946 appears to be the impact 
of a high- velocity HI cloud and/or hyp ernova explo- 
sions (|Elmegreen. Efremov fc LarsenI |2000[ ). The Galaxy 
may harbour similar structures. A possible example is 
the C ygnus superbubble, which contains O B associa- 
tions (jVlemmings. Cordes fc Chatteried |2004 and refer- 
ences therein). 

For a GC, the triggering effects are essentially grav- 
itational. A natural assumption is that GCs, crossing 
the disk every 1 Myr on average, may be responsible for 
some star formation. A possible case relates the origin 
of the GC NGC 6231 to the GC NGC 6397 disk-crossing 
(jRees fc Cudworthl l2003| ). Another possibility is the low- 
mass GC FSR 584 as s tar-formation trigger in the W 3 com- 
plex ()Bica et al.ll2007[ ). 

The OCs Stephenson 2 and B D SB 122 were 

discovere d in 1990 (jStephensonl Il990f) and 

2003 (jBica et al.l I2003D . respectively. 2MASS 
( www. ipac. caltech. edu/2mass/releases/allsky) images 
of both clusters are shown in Fig. [TJ The suspected rich- 
ness of Steph e nson 2 in r ed supergiants was co nfirmed by 
iNakava et all (|200lh and lOrtolani et all (|2002D . providing 
a n age of » 20 Myr, and a distance from the Sun d© = 6 kpc 
fjOrtolani et al. 2002). Both clusters are among the most 
massive OCs known in the Galaxy. Indeed, BDSB 122 has 
14 red supergiants, is located at d© = 5.8 kpc from the 
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Table 1. Present-day cluster positions 



Fig. 1. 5' X 5' 2MASS i^s images of Stephenson 2 (left) and 
BDSB 122 (right). Figure orientation: North to the top and 
East to the left. 
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Fig. 2. The present-day positions (and uncertainties) of 
Stephenson 2, BDSB 122, and w Centauri overplotted on a 
schematic projection of the Galaxy as seen from the North 
pole, with 7.2 kpc as the Sun's distance to the Galactic cen- 
tre. Main structures are identified. 



Sun, an esti mated mass of 2 — 4 x Mq, and an age of 
7 - 12Myr (jFiger et al.ll2006t) . Stephenson 2 has 26 red su- 
pergiants, is located at d0 = 5.8t|!);gkpc from the Sun, has 
an estimated mass o f 4 x 10''' Mq, and an age of 12 — 17Myr 
(jPavies et al.ll2007l) . Their distances from the Sun are the 
same, within uncertainties, and their projected separation 
on the sky is « 1 00 pc. The desig nation Stephenson 2 was 
originally given bv lOrtolani et all 12002 ). and also adopted 
by Dias et al. (2002, and updates). Stephenson 2 and 
BDSB 122 are clearly in the re d supergiant ( RSG) phase 
(|Bica. Santos Jr. fc AllohUllQQOD . [Pavies et all (120071 refer 
to these clusters as RSGC 1 and RSGC2, respectively. 

The positions (and uncertainties) of both clusters, to- 
gether with uj Centauri (NGC 5139), are shown in Fig.[2]su- 
perimposed on a schemati c view of t he Milky Way (based on 
[Momany et al. 2006 and Drimmel fc Sperge l2001.) . The lo- 
cus of BDSB 122 and Stephenson 2 is slightly more internal 
than the Scutum-Crux arm. Several other young clusters 
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from the catalogues of iBica et al.l (|2003l ) and iDutra et al 



(|2003') h ave already b een s tudied in detail (e.g . 'Soare s et al 
[2008; Or tolani et al]|2008t iHanson fc Bubni ck 2008). 

We trace backwards in time the orbits of uj Centauri and 
that of the OC Stephenson 2 (and consequently, also that of 
BDSB 122) in the disk, testing an impact hypothesis for the 
origin of these two massive OCs. Orbit integrations in the 
Galactic potential using as constraints the G C space veloc- 
ity h a ,ve been widely applied to 54 GC s fe.g. iDinescu et all 
120031: lAllen. Moreno fc Pichardoll2008h . 

This paper is structured as follows. In Sect. [5] we study 
the past orbit of Centauri. In Sect. 12.11 the past orbits 
of Stephenson 2 and lo Centauri are compared looking for 
spatial and time coincidence. Conclusions are in Sect. [31 



2. uj Centauri as a projectile 

to Centauri, the m ost massive Galactic GC (4 x 10^ Mq- 
iNakava et al.ll200l}) . has a metallicity spread and a flat den- 
sity distributio n typical of a dwarf gala xy nucleus captured 
by the Galaxy (jBekki fc Freemanll2003f ). Thus, irrespective 
of the existence of young massive clusters, somehow associ- 
ated with the impact site, the orbit of u Centauri under the 
Galactic potential is worth studying, especially to look for 
consequences of the last disk passage. Evidence of a similar 
disk impact and a star-forming event has been observed in 
the spiral galaxy NGC 4559 with Hubble Space Telescope 
(HST), XMM-Newton, and ground-based (SCP05) data. 
The age of the star-forming complex is ^ 30 Myr, with a 
ring-like distribution. It appears to be an expanding wave 
of star formation, triggered by an initial density pertur- 
bation. The most likely triggering mechanism was a col- 
lision with a satellite dwarf galaxy crossing through the 
gas-rich outer disk of NGC 4559, which appears to be 
the dwarf galaxy visible a few arcsec NW of the com- 
plex. The picture is remi niscent of a scaled-down version 
of the Cartwheel galaxy jStruck-Marcell fc HigdonI 119931 : 
IStruck-Marcell et al.|[T996l ). 

As another example, proper motions (PMs) and ra- 
dial velocity suggest that the GC NGC 6397 crossed 
the Galactic disk 5 Myr ago, which possibly triggered 
the formation of the OC NGC 6231 (|Rees fc CudworthI 
200 ^, thus lending support to the present scenario 
( W allin. Higdon fc Stavelev-SmithI [l99l . NGC 6397 and 
NGC 6231 are closely projected on the sky {A£ 5°, 
A6 ~ 13°). However, in the case of w Centauri as possi- 
ble generator of BDSB 122 and Stephenson 2, the GC is 
now widely apart from the pair of massive OCs (A^ w 77°, 
A6 w 15°). Thus, PM and radial velocity are fundamental 
constraints for the analysis of lo Centauri, and impact solu- 
tions require a detailed integration of the orbit across the 
Galactic potential. 
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Fig. 3. Top-left panel: Galactocentric XY-plane projection of the Centauri orbit over the past 2Gyr. Top-right: The 
past 30Myr orbit of wCentauri (solid line) for i?GC ~ 7.2 kpc. Additional neighbouring orbits (dotted) are, from bottom 
to top: JHB96 (-10%), JHB96 (-1-10%), FSC96. The impact site on the disk is shown by the empty square. Arrows 
indicate orbit direction. Orbits of Stephenson 2 for the assumed distance from the Sun (and uncertainties) are shown. 
The corresponding XZ and YZ projections are in the bottom panels. Empty symbols over the Stephenson 2 orbits indicate 
its possible positions 24Myr ago. The Sun at its present position (asterisk) and the Galactic Centre (filled circle) are 
shown. 



2.1. Orbit computation 

The present three-component mass-distribution model 
of the Galaxy follows that employed in the study 
of a high-velocity black hole on a Galactic-ha l o or- 
bit in the solar neighbourhood (|Mirabel et al.l l200ll 
and references therein). In short, we use the three- 
compo nent model of iJohnston. Hernquist fc Bolt3 
(|l996t) - hereafter JHB96 - in which the disk, 
spheroidal, and halo gravitational potentials are de- 
scribed as 0disk(i?, z) = -GMdisk/\J R"^ + {a + ^/ z"^ + h'^f 
(iMivamoto fc Naeail [Hzl), (^sphcr(i?) = -^^Sf^ 



(|Hernauisti fT99(Th. and </.haio(i?) = vl^^^ln {R^ + cP) , 
where Mjisk = 1 x 10" Mq, M.phor = 3.4 x lO^o Mq, 
'^haio = 128kms~^, R and z are the cylindrical coordi- 
nates, and the scale lengths a = 6.5 kpc, h = 0.26 kpc, 
c = 0.7 kpc, and d = 12.0 kpc. Table [1] shows the Galactic 
and Equatorial coordinates of t he thr ee clusters. Following 
iMizutani. Chiba fc Sakamotcl (|2003[ ). the relevant pa- 
rameters for computing the motion of lo Centauri are 
the distance from the Sun ~ 5.3 ± 0.5 kpc, the PM 
components (masyr~^) /Uq cos((5) = — 5.08 ± 0.35 and 
lis = —3.57 ± 0.34, and finally the heliocentric radial 
velocity Vr = 232.5 ± 0.7 km s'^ 
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Fig. 4. Close-up of the impact site. Left panel: Orbits computed with Rgc = 7.2 kpc (empty symbols) and Rgc — 7.6 kpc 
(filled symbols). Right: Same as left panel but including error distribution for the lu Centauri impact site and Stephenson 2 
proto-cluster position. A random selection of impact sites (open circles) is shown within the uj Centauri error ellipsoid. 



The models were computed with Rqc = 7.2 kpc 
(jBica et al.ll2006l ) as the distance of the Sun to the Galactic 
centre. The Galactic velocities of u) Centauri are U = 
54.3 ± 9.5kms-i, V = -44.2 ± 8.2 kms'^, and W = 
— 1.3 ± 13.0km s^^. Alternatively, we al so computed or- 
bits w ith i?GC = 7.6 kpc, obtained bv lEisenhauer et "all 
(|2005( ). It should be noted that recently, by means of sta- 
tistical parallax of central sta rs, Trippe et al . (2008) found 
Rgc = 8.07±0.32kpc, while iGhez et all (|2008h . with the 
orbit of one star close to the black hole, found Rgc = 
8.0 ± 0.6 kpc or Rgc — 8.4 ± 0.4 kpc, under different as- 
sumptions. Cluster distances are heliocentric, which do not 
depend on Rgc', on the other hand, the value of Rgc has 
some effect on the potentials, which thus, affects orbit com- 
putation. Since the difference between the adopted value of 
Rgc and the more recent ones is not exceeding, the value 
of Rgc should not influence much the present results. 

Based on the rotation curves of lBrand fc Bhti (|l993l) 
and Russcil (2003), and an estimate with the galaxy mass 
model described above ([Mirabel et al.l [200l[ ). we derived 
Vc = 214 ±4 km s~^ as the orbital velocity of Stephenson 2. 
The nearly flat Galactic rotation curve at the Stephenson 2 
position allows to adopt this circular velocity also for the or- 
bits corresponding to distance uncertainties (Sect. [T]). The 
orbit of IV Centauri, c omputed back over 2 Gyr, is co i npara - 
ble to that derived bv lMizutani. Chiba fc Sakamotol (|2003t) . 
in particular the Rosette pattern seen projected on the XY 
plane (Fig. [3]). The simulation shows that Centauri hit 
the disk as recently as 24 ± 2 Myr ago. Because of such a 
short time, fossil remains of this event may be detectable 
in the disk nowadays. 

Figure [3] (top-left panel) shows the Galactic XY-plane 
projection of the last 2 Gyr orbital motion of co Centauri. 
In the remaining panels we focus on the last 30 Myr of the 
motion of w Centauri and Stephenson 2. For Stephenson 2 
we consider the different orbits that result from the 
adopted distance from the Sun and corresponding uncer- 



tainties (Sect. [T|). The interesting fact is that the orbit of 
Stephenson 2 passes close to the impact site of lu Centauri at 
a comparable time, within uncertainties (see below). Since 
Stephenson 2 and BDSB 122 have almost the same position 
(within uncertainties), the same conclusions hold for the 
latter cluster. The XZ and YZ-plane projections (bottom 
panels) show that uj Centauri emerged at « 45° from the 
plane to its present position. 

To probe orbital uncertainties owing to the adopted 
potential, we also employed the potenti al model of 
iFlynn. Sommer-Larsen fc ChristensenI (|l996l ) - hereafter 
FSC96 - and tested consequences of variations of ±10% in 
the input parameters of JIIB96. The results are shown in 
Fig. [3] (top-right panel), from which we conclude that orbit 
variations due to the adopted potential are much smaller 
than our error ellipsoid (Fig. 31 right panel). 

Close-ups of the uj Centauri impact site and the back- 
traced positions of Stephenson 2 are shown in Fig. 2] (left 
panel) for a Sun's distance to the Galactic centre of 7.2 kpc 
and 7.6 kpc. It is clear that the latter value is not critical 
for the encounter. The right panel shows the error ellip- 
soid of several impact site simulations computed by vary- 
ing initial conditions according to the errors in the differ- 
ent relevant input quantities. The ellipsoid reflects varia- 
tions implied by velocity uncertainties in the PM, radial 
velocity and present position of iv Centauri along the line of 
sight in the (U,V,W) velocities. The impact obtained with 
a Galactocentric distance Rgc = 7.6 kpc is also shown. 
The disk-orbit of Stephenson 2 crosses the uj Centauri ellip- 
soid error distribution. The range in impact site to proto- 
cluster separations contains distances smaller than « 1 kpc, 
with an average separation of ~ 500 pc (intersection area in 
Fig. [31 right panel) . Larger separations would require pro- 
hibitive expansion velocities, despite the fact that we are 
dealing with an encounter in a denser, central part of the 
disk, while in NGC4559, the event was external. 
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For the GC-induced formation hypothesis to be valid, 
the time-scales associated with the onset of star for- 
mation (after impact), duration of star formation and 
the cluster age, shoul d be compatible with the disk- 
crossing age. Following [Putte fc Cropped (|2009D . the first 
time scale in not well known, ranging from virtually 
instantaneous, i.e. negligible as compared to the clus- 
ter age, to 15Myr CMpine_fc_Du yert _199j) and 30Myr 
(|Wallin. Higdon fc Stavelev-Smithlll996h . The star forma- 
tion time-scale may b e short, « 2 x 10^ yr, as suggested by 
iMcKee fc Tall (|2002D for stars more massive than 8 Mq . 
Given that the ages of Stephenson 2 and BDSB 122 are 
within 12-17 Myr and 7-12 Myr, respectively, w Centauri, 
which crossed the disk « 24 Myr ago, may have triggered 
their formation only if the star-formation onset occurred in 
less than w 15 Myr, which is within the accepted range. In 
the case of NG C4559 these time -scales combined are less 
than ~ 30 Myr (jSoria et al.ll2005D . 

The above clues suggest an interesting event, that the 
most recent crossing of co Centauri through the disk oc- 
curred very close to the sites where two massive OCs were 
formed. Both Stephenson 2 and BDSB 122 are somewhat 
younger than the age of the impact, and the differences of 
a few Myrs are consistent with the shock propagation and 
subsequent star formation. The overall evidence gathered 
in the present analysis supports uj Centauri as the origin of 
this localised star formation in the Galaxy, which harbours 
two of the more massive known OCs. 

This work suggests a scenario where the disk passage 
of GCs can generate OCs, massive ones in particular, as 
indicated by the orbit of w Centauri and its impact site. 
As a consequence, not all OC formation is induced by the 
spiral density wave mechanism in spiral arms. 

3. Summary and conclusions 

Globular clusters orbiting the bulge and halo of the Galaxy 
cross the disk on average once every 1 Myr, and significant 
physical effects of such events on the disk are expected to 
occur. For instance, the impact of a GC passing through 
the disk can trigger star formation either by accumulation 
of gas clouds around the impact site or by the production of 
an expanding mechanical wave. Time delays are expected in 
both cases because of collapse and fragmentation of molecu- 
lar clouds before star formation. This phenomenon has been 
recently observed in the galaxy NGC4559 (e.g. SCP05). If 
this mechanism operates frequently in the Galaxy, the most 
massive GC u Centauri can be taken as an ideal projectile 
to prospect the state of its last impact site in the disk. 

As w Centauri hit the Galactic disk « 24 Myr ago, a 
major star-forming event appears to have occurred close 
(^ 1 kpc) to the locus that generated two of the most 
massive young OCs known in the Galaxy, BDSB 122 and 
Stephenson 2. We suggest the connection of these events 
in a way similar to the impact and shock wave observed 
in NGC4559 (SCP05). We use a model of the Galactic 
potential to integrate the orbit of to Centauri. As shown 
in Fig. 31 when the uncertainties in space velocity, dis- 
tances, and potential are considered, the error-distributions 
of oj Centauri impact site and the birth-site of Stephenson 2 
overlap. Such overlap suggests a scenario where the disk 
passage and formation of the pair of OCs may be phys- 
ically connected. Alternatively, the time coincidence may 
have occurred within a separation ^ 1 kpc. In such case. 



the expanding bubble scenario such as that in NGC 4559 
would apply. The latter case is more probable, since two 
clu sters have be en formed. 

iLevvl (|2000t ) performed 2D hydrodynamic simulations 
to study the impact of GCs on the Galactic disk in the 
presence of available gas. They found that the moving GC 
causes a shock wave in the gas that propagates through the 
disk in a kpc scale, thus produc ing star formation. More 
recently, iPutte fc Cropper] (|2009D simulated in detail the 
effects of GC impacts on the disk that basi cally c onfirm 
IWalhn. Higdon fc Stavelev-SmithI (119961) and lLevvl (|2000D . 
even in the absence of gas at the impact site. They found 
focussing of disk material with compression on a scale of 
~ 10 pc, that may subsequently attract gas leading to star 
formation. The compression increases with the GC mass. 

At this point an interesting question arises. For a rate 
of ~ 1 GC impact per Myr, it can be expected a high prob- 
ability of occurring one GC impact within 1-2 kpc of any 
location within the inner Galaxy in about 10 Myr. However, 
the star-formation efficie ncy of thes e even ts appears to be 
low, according to Putte fc Cropp"e^ (|2009[ ). Indeed, of the 
54 GCs with accurate proper motions studied by them, only 
three appear to be associated to young OCs. Possibly, con- 
ditions like GC mass and impact site properties, such as 
the availability of molecular gas, temperature and density, 
constrain the star-formation efficiency. 

Evidence drawn from the present work suggests that 
GCs can be progenitors of massive OCs. In particular we fo- 
cused on UJ Centauri. Density-wave shocks are possibly not 
the sole responsible for the formation of the more massive 
OCs, a possibility to be further explored, both theoretically 
and observationally. 
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